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Saturation in High-energy QCD 

Gregory Soyez* 

SPhT, CEA/Saclay, Orme des Merisiers, Bat 774, F-91191 Gif-sur-Yvette cedex, Franc^ 

In these proceedings, I shall review the basic concepts of perturbative QCD in its high-energy limit, empha- 
sising the approach to the unitarity limit, usually referred to as saturation. I shall explain the basic framework 
showing the need for saturation, first, from a simple picture of the high-energy behaviour, then, giving a short 
derivation of the equation driving this evolution. In the second part, I shall exhibit an analogy with statistical 
physics and show how this allows to derive geometric scaling in QCD with saturation. I shall finally consider 
the effects of gluon-number fluctuations on this picture. 



I. INTRODUCTION 

The quest for the high-energy behaviour of perturbative 
QCD started thirty years ago, soon after QCD was proposed 
as the fundamental theory of strong interactions. To grasp this 
problem, it has been realised that saturation effects were to be 
considered in order to satisfy unitarity constraints. Through- 
out these proceedings, I shall summarised the present status of 
our understanding of this limit, which basically means finding 
an equation giving the evolution of amplitudes towards high 
energy and finding the properties of its solutions. 




log(Q 2 ) 

FIG. 1: Picture of the proton in DIS 

In order to motivate the physical picture we want to repro- 
duce, let us consider the problem of Deep Inelastic Scattering 
(DIS) i.e. Y*P ~ ¥ X collisions. This process depends on two 
kinematical variables: the virtuality Q 2 of the photon, mea- 
suring the resolution at which the photon scans the proton, 
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and the Bjorken x, given by the fraction of the proton mo- 
mentum carried by the parton struck by the virtual photon in 
a frame where the proton is moving fast. This last variable is 
related to the centre-of-mass energy s = Q 2 /x, meaning that 
the high-energy limit corresponds to the small-x limit. For 
practical purposes, we introduce the rapidity defined through 

y = io g (i/x). 



In Figure ^ we have represented the typical configuration 
of the proton in different domains of this phase space. We 
start with a proton at low Q 2 and low energy, represented as 
three partons (bottom-left part of the picture). If we increase 
the virtuality of the photon, we are able to resolve its partons 
into smaller ones. This type of evolution is described in QCD 
by the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) 
equation []]]. In that limit, the number of partons rises log- 
arithmically while their typical size decreases like l/Q 2 so 
that the proton becomes more and more dilute. In what fol- 
lows, we shall be interested in the other type of evolution, 
namely the evolution towards smaller values of x for constant 
Q 2 . If one boosts the proton, what typically happens is that 
we create new partons of a size comparable to the parents' 
ones. Obviously, during this evolution, known as the Balitsky- 
Fadin-Kuraev-Lipatov (BFKL) evolution |2], the proton tends 
to the black-disc limit for which partons start to overlap. At 
that level, they start to interact among themselves and unitar- 
ity corrections are to be considered 1 3]. This blackening of the 
proton is called saturation and is what we shall describe from 
perturbative QCD in these proceedings. 



We shall begin our analysis with a description of the BFKL 
equation. This naturally leads to the introduction of unitar- 
ity corrections and to new evolution equations - the Balit- 
sky/JIMWLK hierarchy jj, |5|] and the Balitsky-Kovchegov 
(BK) equation |4, 6]. We shall then show that, within relevant 
approximations, the BK equation reduces 0] to the Fisher- 
Kolmogorov-Petrovsky-Piscounov(F-KPP) equation |8|, well 
studied in statistical physics. This allows for a description 0] 
of the asymptotic behaviour of the scattering amplitudes in 
terms of travelling waves, translating in QCD into geomet- 
ric scaling In the last section, we shall discuss the ef- 
fect of g luon-number fluctuations, which have recently proven 
JToL ll ll IT2I fl3ll to have important consequences on the ap- 
proach to saturation. 
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II. EVOLUTION TOWARDS HIGH ENERGY 

Let us start with a fast-moving quark of momentum p. This 
quark emits Bremsstrahlung gluons characterised by their 
transverse momentum k 2 and their longitudinal momentum 
k z — xp. The probability for emitting such a gluon is 



dP> 



dx dk 2 , 
x k^ 



When we consider large transverse momentum (high Q in 
DIS), the collinear divergences dk\ /k 2 need to be resummed 
and this gives rise to a renormalisation group equation known 
as the DGLAP equation. Throughout this proceedings, we 
shall instead consider the limit of high energy, in which, we 
are sensitive to the gluons of small longitudinal momentum 
j « 1. This leaves a large phase-space to have successive 
gluon emissions satisfying 1 ^> x\ 3> ■ ■ ■ S> x„ 3> x which 
gives a contribution of order 



a: 



dx\ 



1 



a"log(l/x)". 



Even when the coupling a s is small enough to ensure appli- 
cability of perturbation theory, for small values of x i.e. high- 
energy, a s log(l fx") ~ 1 and all these contributions have to be 
resummed. 



FIG. 2: Gluon emission in the dipole picture 

In order to simplify the discussion, let us consider the case 
of a large number of colours. In that limit, a gluon of trans- 
verse coordinate z can be considered as a quark- an tiquark pair 
at point z. In this framework, it is convenient to consider the 
evolution of a colourless quark-antiquark dipole 0. Indeed, 
at high-energy, all possible gluon emissions from one dipole 
are tantamount to this dipole splitting into two child dipoles 
as depicted in figure[2] We can thus forget about gluons them- 
selves and talk only in term of dipoles. 

Let us denotes by (T xy ) the scattering amplitude 1 for a 
dipole made of a quark at transverse coordinate x and an an- 
tiquark at y. Getting the evolution of (T) from the dipole pic- 
ture is pretty straightforward. Indeed, if one boosts the dipole 
(x,y) from a rapidity Y to a rapidity Y + 8Y, this dipole shall 
split into two dipoles (x,z) and (z,y). Each of those dipoles 
can interact with the target which leads to the following evo- 
lution equation 2 for (T) 



a 

(T xy ) = — 
V xy/ 271 



M xyz ((T xi 



(1) 



: : 
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o o 














FIG. 3: Multiple scattering between a target and a projectile made 
of two dipoles. This can be interpreted as pomeron merging in the 
target. 



with a = a s N c /n and 



5W xyz = 



(x-z) 2 

(x-y) 



(y-z) 2 



(x-z) 2 (z-y) 2 

is the probability density for dipole splitting 3 . Equation Q is 
the BFKL equation, derived in the mid-seventies |2] by Balit- 
sky, Fadin, Kuraev and Lipatov. 

This equation being linear, one expects its solution to 
grow exponentially. To be more explicit, if one neglects 



impact-parameter dependence, i.e. assumes that (J, 
(T(r = |x — y|)), the solution of Q is 



(T(r)} 



dy 

2/7C 



ToO) 



aW-jlogirl/r 2 ) 



(2) 



where = 2\|/(1) — \|/(y) — v|/(l - y) and Tq(j) describes 
the initial condition. For a fixed dipole size, one finds the high 
energy behaviour in the saddle point approximation 



(T(r)) -exp 



(tipY- 



2x"(i/2)aF 



with <Mp = 47tlog(2)a being the BFKL pomeron intercept. 
This expression suffers from two major problems: first, even if 
one start with an initial condition peaked around a small value 
of r, when rapidity increases, the amplitude starts to diffuse to 
the large dipole sizes i.e. to the non-perturbative domain. Sec- 
ond, the solution of the BFKL equation grows exponentially 
with rapidity. It hence violates the unitarity constraint T < 1 
obtained from first principles. 

This unitarity problem is precisely the point where we meet 
the requirement to take into account multiple interactions . The 
most straightforward way to see this is to come back to the 
splitting of a dipole (x,y) into (x,z) and (z,y). When the tar- 
get becomes dense enough, both dipoles can scatter on it. We 
thus have to take into account this contribution, represented 
in figure [5] This leads to a quadratic suppression term in the 
evolution equation which becomes 

h (Txy) = ^ «r„) + <r zy ) - (T xy ) - (r x 2 z;zy )) . (3) 



1 The notation {■) denotes the average over all possible realisations of the 
target colour field. It can be understood as the expectation value for the 
T-matrix operator over the target wavefunction. 

2 The last term, coming with a minus sign, corresponds to virtual corrections 



3 The first line shows explicitly the two contributions coming from emission 
by the quark and the antiquark as depicted in fieureEl 
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The first thing to remark is that this new term is of the same 
order as the previous ones when T 2 or, equivalently, when 
T ~ 1. It is thus, as expected, a mandatory contribution near 
the unitarity limit. However, equation Q involves a new ob- 
ject, namely (r xz - zy ) which probes correlations inside of the 
target. In a general framework, one should then write down an 
equation for {T 2 ), which will involve {T 2 ) through BFKL- 
like contributions and (r 3 ) from unitarity requirements. This 
ends up with a complete hierarchy, giving the evolution for 
each (T k ), known as the (large-jV c ) Balitsky hierarchy @|. 

If the target is sufficiently large and homogeneous, one can 
simply assume {T 2 z . zy ) = (T xz ) (T zy ), ending up with a closed 
equation for (T) 



W Y [p] 



W Y [p] 



XV/ 



a 
2% 



^xyz ((7x Z ) + (T zy ) - (T xy ) - (T xz ) (T zy )) 



(4) 

This last expression is the Balitsky-Kovchegov (BK) equa- 
tion |H,|(|, which is the most simple equation one can obtain 
from perturbative QCD by including both the BFKL contribu- 
tions at high-energy and the corrections from unitarity. It is 
easy to check that (T) = is an unstable fix point of the BK 
equation, while (T) = 1 is a stable fix point, ensuring unitar- 
ity is satisfied. In general, unitarity effects cut the increase of 
the amplitude for scales smaller than a typical scale, called the 
saturation momentum, increasing with rapidity. In particular, 
it solves the problem of diffusion to the infrared obtained with 
linear BFKL equation. We shall come back with more details 
to the study of the solution of the BK equation for the next 
section. 

If one wants to study those effects beyond the large-A^ £ 
approximation, one better used the colour-glass-condensate 
(CGC) formalism. This is an effective theory in which one 
consider the radiation of soft gluons in a strong background 
field which introduces the nonlinear effects requested at satu- 
ration. At a given rapidity, the target is considered as a random 
colour source p u (x) and radiates a field according to the QCD 
Yang-Mills equations. We then study the rapidity evolution of 
the functional Wy [p] giving the probability to find a source p 
at rapidity Y. This functional allow to compute the average 
value of an operator using 



2>[p]Wy[ P ]0[p]. 



(5) 



We thus need to find the evolution of W [p] . One step of 
evolution is represented diagrammatically in figure 0] As for 
the BFKL equation, one soft gluon is radiated (yielding a 
oc v log(l/jc) factor). The difference comes from the fact that 
when the target is sufficiently dense, one has to dress the prop- 
agator of this soft gluon with the field created by the target. 
These emissions are then reabsorbed inside the probability 
W[p] in order to define it at rapidity Y + 8Y. Each background 
gluon insertion gives a factor gu with g the QCD coupling 
and a the A + colour field of the target. Including these ef- 
fects becomes mandatory when the field is of order 1/g. The 
amplitude is defined through 



(7xy) = l-^(tr(w x t W y ) 





FIG. 4: One step of JIMWLK evolution: the emitted soft gluon prop- 
agating in the string target field is used to redefine the probability 
W[p] at rapidity Y + 8Y. The case of two gluons only corresponds to 
BFKL. 



where W x is the Wilson line operator 



W x = P exp 



dx oc(.x ,x) 



which means that in the leading order, (r xy ) °< g (oc x — 0C y ) . 
This shows that the strong field a ~ 1 jg corresponds to T ~ 1 
and thus to saturation corrections. A careful treatment of the 
diagrams in figureheads to the JIMWLK equation 



oywv[p] = 



(7) 



(6) 



with the kernel % computable in perturbative QCD. 

This functional integro-differential equation gives the ra- 
pidity evolution of W [p] . Using relations like Q or 0, one 
can obtain the evolution of the scattering amplitudes. The 
resulting equations are equivalent to the Balitsky hierarchy 
yj. These equations going beyond large N c , they involves 
non-dipolar operators. For example, the equation for two 
dipoles depends on a sextupole operator (build out of six Wil- 
son lines). The formalism beyond large N c being rather tricky, 
in what follows we shall stick to the BK equation, which is 
sufficient to grasp the basic concepts of saturation. 

Very recently, it has been realised that the question of the 
high-energy limit of QCD, even at the leading logarithmic ac- 
curacy, is not fully described by the Balitsky hierarchy. To 
see this more precisely, let us consider in more details the 
evolution equation for {T 2 ). It corresponds to the scatter- 
ing of two dipoles off a target and the diagrams contributing 
to one step of rapidity evolution are shown in figure [5] The 
first graph on the left of this figure is simply the BFKL-ladder 
contribution for which dy (T 2 ^ °< (T 2 )- The diagram in the 
middle corresponds to multiple scattering (see also figure^. 
It accounts for unitarity corrections and gives a contribution 
dr (T 2 ) x (T 3 )- If those two graphs were the only relevant 
ones in the evolution, one would obtain the (large-A^) Bal- 
itsky equation. Nonetheless, it has been argued that a third 
contribution, represented by the rightmost diagram of figure 
[5] has to be included. This graph describes multiple scattering 
off the projectile or, equivalently, gluon-number fluctuations 
inside the target. After the pomeron-merging term, contribut- 
ing to saturation, one thus now also includes the pomeron- 
splitting contribution, which gives rise to a new term in the 
evolution equation 3y(r 2 ) °< (T). For large-A^ £ . and at the 
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FIG. 5: Diagrams contributing to the evolution of (T 2 ). They correspond, from left to right, to BFKL ladders, saturation corrections and 
fluctuations effects. 



two-gluon-exchange level, this new term has been computed 
O and leads to 



1 a /a s \ 2 



(8) 



iyi;x 2 y 2 /ifluct 2 2n\2n 

7uv) ~t~ (1 2), 



with Aq being the dipole-dipole amplitude at the two-gluon- 
exchange level 



^o(xy|uv) = - log 



(x-u) 2 (y-v) 2 
(x-v) 2 (y-u) 2 



This object naturally comes out when we relate the scattering 
amplitude with the dipole density «: 



a; / J3 (xy|uv) («„ 



(9) 



This relation can be inverted to obtain (n 



XV'/ 



The correction (|8} becomes of the same order as the BFKL 
contribution when T 2 ~ a 2 7\ or T ~ a 2 , i.e. in the dilute 
regime where, as expected, fluctuations should lead to im- 
portant effects. Although this may at first sight seems irrel- 
evant for the physics of saturation, one has to realise that, 
due to colour transparency, the amplitude becomes arbitrar- 
ily small for small dipole size, hence, even for dense ob- 
jects, there is always some region of the phase-space where 
the system is dilute and its evolution governed by fluctua- 
tion effects. Once the system has evolved through a fluctu- 
ation, it grows, from BFKL evolution, like two pomerons, 
i.e. T ~ a 2 exp(2(0/>5 / ). This has to be compared with the 
one-pomeron exchange T ~ exp((x)pY). We see that at very 
large energies Y > co^ 1 log(l/oc 2 ), BFKL evolution compen- 
sates the extra factor of oc 2 coming from the initial fluctuation. 
We shall see more precisely in the last section that those fluc- 
tuation effects in the dilute tail have important consequence 
on the saturation physics. 

Finally, let us conclude this section by quoting that the 
equation (|8j, including linear, saturation and fluctuation ef- 
fects is the most complete equation known in perturbative 
high-energy QCD so far. Its extension beyond the large-A^- 
approximation is still a challenging problem. 



III. PROPERTIES OF THE SCATTERING AMPLITUDES 

In this section, we shall sketch out the main properties of 
the solutions of the evolution equations towards high energy. 
As we shall explain in detail in the next lines the major source 
of information comes from the analogy between the BK 
equation and the Fisher-Kolmogorov-Petrovsky-Piscounov 
(F-KPP) equation which is very well studied in statistical 
physics. This section shall be devoted to derive this anal- 
ogy and its consequences in the case of the BK equation. We 
leave the discussion concerning the effect of fluctuations for 
the next section. 

So, let us start with the BK equation (0). We shall first re- 
strict ourselves to the impact-parameter-independent version 
of the equation, for which one can easily go to momentum 
space by using 4 



T(k) 



1 

271 



d 2 r 



ik.r 



T(r) 



dr 



J (kr)T(r) 



(10) 



and where the equation becomes 

d Y T(k) = ax(-d L )T(k)-aT 2 (k), 

with L = log^/k 2 )), ko being a soft reference scale. 

In order to simplify the problem, we shall work in the sad- 
dle point approximation (often referred to as the diffusive ap- 
proximation) and expand the BFKL kernel to second order 
around y = 1/2, transforming the complicated differential op- 
erator in equation (II 01 into a second-order operator 



x(-d L )T(k)n % a)+ix"(k)(d L 



Up to a linear change of variable switching from Y and L to 
time t = uY and space x = L + cst.F and a renormalisation of 
the amplitude T — > u = cst.T, the BK equation in the saddle 
point approximation becomes 



d t u(x,t) = d 2 u(x,t) + u(x,t) 



■ u 2 (x,t) 



This equation is nothing but the F-KPP equation studied in 
statistical physics since forty years and applying to many 



4 In this section, we omit the (•) quotations as they do not play any role in 
the mean-field approximation. 
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lowing conditions are satisfied 

1 . the amplitude has as unstable fix point and 1 as stable 
fix point, 

2. the initial condition decreases faster than exp(— y c L) at 
large L (see below for the general definition of y c ), 

3. the equation obtained by neglecting the nonlinear terms 
admits superposition of waves as solution. 

Then, travelling waves are formed during evolution with 
the critical parameters obtained, from the linear kernel only, 
through the relation 



FIG. 6: Numerical simulation of the rapidity-evolution of the BK 
equation. We start at Y = from the leftmost amplitude and evolve 
to higher Y using fTol . Amplitude is shown for Y = 5, 10, 15,20,25 
and clearly exhibits a travelling-wave pattern. 



So 
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FIG. 7: Rapidity evolution of the saturation scale extracted from fig- 
ure!^] More precisely, the speed of the wave, i.e. the exponent of the 
saturation scale, is plotted. It goes to a constant as it should. 



problem. It describes reaction diffusion processes in the 
mean-field approximation where one can have creation and 
annihilation of particles locally and diffusion to neighbouring 
site e.g. the density of population in the centre of Rio. 

One knows that the F-KPP equation admits travelling waves 
as asymptotic solutions i.e. there exists a critical velocity v c 
and a critical slope y c , determined only from the knowledge 
of the linear kernel of the F-KPP equation d 2 . + 1, such that 
(see figure [5] for a pictorial representation coming from the 
BK equation which we shall explain a bit later) 



u{x,t) 



u(x - v c t) e ->(*-^) log(x - Vc t). 



Although working in the diffusive approximation makes the 
discussion easier due to the direct analogy with the F-KPP 
equation, this assumption is not required. Indeed, for an equa- 
tion with a more complicated kernel (than d 2 + 1), one can 
prove the existence of travelling waves provided the three fol- 



One can show that this corresponds to the selection of the 
wave with the minimal speed x(Y)/Y and the point where 
group and phase velocities are equal. 

In the case of the BK equation, the first condition is satisfied 
due to BFKL growth and saturation, the second comes from 
colour transparency (T « 1 jk 2 = exp(— L)) and the last one is 
equivalent to the solution 10 (with r/ro replaced by ko/k) of 
the BFKL equation obtained by neglecting nonlinear terms in 
EJ. 

The critical parameters (for the complete BFKL kernel i.e. 
for (|T0|l), are found to be y c « 0.6275 and v c « 4.8836a. Writ- 
ten in terms of the QCD variables Y and k 2 , one can then write 
the asymptotic solution for the impact-parameter-independent 
BK equation under the form 



T(k;Y) 



-Yc 



log 



QW) 



with the saturation scale given by 



Q 2 (Y) r =~k 2 exp 



v,y--log(Y) 



(11) 



(12) 



The formation of a travelling-wave pattern when energy in- 
creases can easily be seen on numerical simulations of the BK 
equation. As displayed in figure |6j if one start with a steep 
enough initial condition (leftmost curve), the amplitude in- 
creases with energy and a wave moving into the dilute domain 
gets formed. From that simulation, one can extract the satura- 
tion scale by solving T(Q 2 , Y) = N at each values of Y and for 
a fixed threshold N. In figure[71 we have plotted \og{Q^)/Y 
which goes to a constant value as it is expected from d!2i . 

Equation dl It has a remarkable property: it proves that at 
high energy the amplitude, a priori a function of both Y and 
k, depends only on the ratio between k and the saturation mo- 
mentum. This property, known as geometric scaling, has been 
observed 01 in the HERA measurements of the proton struc- 
ture function. The fact that geometric scaling can be derived 
from the BK equation is one of the most important indication 
for the experimental observation of saturation. The saturation 
scale obtained from the structure function data is of the order 
of 1 GeV for x ~ 10~ 5 . This means that high-energy QCD can 
indeed be studied from perturbation theory. 
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Up to now, we have only discussed the impact-parameter 
independent BK equation. We might therefore ask whether 
or not these arguments extend to the full equation including 
all phase-space degrees of freedom. A dipole of transverse 
coordinates (x, y) is then represented through its size r = x y 
and impact parameter b = (x + y)/2. The problem is then that 
dipole splitting is non-local in impact parameter. This means 
that the BK equation couples different values of b and we can 
not apply directly the previous arguments for each value of b. 
Again, the solution consists in moving to momentum space 
and replace the impact parameter by momentum transfer q 



j2„ j2„ „ik.x ;(q-k).y Tfcj) 



T (k,q) = / d 2 xd 2 ye ,kx e ,((i - k > y 



(x-y) 2 



with 



A 2 = i*o if *»> » 9 
1 q 2 if q > ko 

It is interesting to notice that the critical slope y c ad speed v e , 
obtained from the BFKL kernel only, are the same as for the 
impact-parameter-independent case. 

The important point at this stage is that this study predicts 
geometric scaling at non-zero momentum transfer. Experi- 
mental measurements of the DVCS cross-sections or diffrac- 
tive p-meson production are good candidates and can provide 
more evidence for saturation. 



The BK equation then takes a form 1 15] for which the BFKL 
kernel is local in q 



3rf(k,q) 

a r d 2 ]<! 

a 

27t 



(13) 



(k-k'f 



?(k',q) 



k' 2 



(q-k) 



2 I 



(q-k') 2 



J rf 2 ff(k,k')f(k-k',q-k') 



We can now proceed with this equation in a similar way as 
for the the impact-parameter-independent case. The existence 
of travelling waves required the three conditions stated previ- 
ously to be satisfied. The two first ones are trivially satisfied 5 . 
For the third condition to be valid, we need to find solutions 
of the linear part of (II 31 , i.e. the complete BFKL equation, 
which can be expressed as a superposition of waves. This is 
merely technical as it consists in a careful treatment of the so- 
lutions of the BFKL equation 1 16] in its full glory, so I shall 
directly jump to the results. We can show 11711 that when the 
dipole momentum k is much larger than the momentum trans- 
fer q and the typical scale of the target ko, the solutions of the 
BFKL equation are a superposition of waves and therefore, 
we obtain travelling waves for the full BK equation. More 
precisely, the high-energy behaviour of the amplitude takes 
the form 



T(k,q;Y) 



Q 2 (q;Y) 



log 



Qi{q;Y) 



This expression is exactly the same as for the previous case 
except that the saturation scale now depends on momentum 
transfer 



g^T) = A 2 exp 



v c y--io g( y) 



5 Indeed, we have saturation and BFKL growth ensuring the first condition 
and colour transparency still apply for the second one. 



IV. HIGH-ENERGY QCD WITH FLUCTUATIONS 

Let us now consider the effect of adding the fluctuation con- 
tribution (|8}. This new term, added to the full Balitsky hier- 
archy, turns a simple non-linear equation into a infinite hier- 
archy with a complicated transverse-plane dependence 6 . We 
can however simplify the equation to a tractable problem by 
performing a local-fluctuations approximation. This amounts 
to simplify the dipole-dipole scattering amplitude j?q used to 
relate the dipole density n with the scattering amplitude T. 
We assume that two dipoles interact if they are of the same 
size and if their centre-of-mass are sufficiently close to allow 
for an overlap of the two dipoles. Within this approximation 
equation (|9} is replaced 



(r xy )=Ka2|x-y|- 



'xy/ > 



where K is an unknown fudge factor of order 1 . 

Once this is done, the remaining steps are no particularly 
illuminating: one Fourier- transform the dipole size r = |x — y | 
into momentum k and perform a coarse-graining approxima- 
tion to get rid of the impact-parameter dependence of the fluc- 
tuation term. This computation results into the following sim- 
plified form of the hierarchy (we give only the expression for 
(T) and (r 2 ) for simplicity) 



dr (T k ) 



= ax(-d L )(T k )-a(T 2 k ), 

= a X (-d Ll )(T 2 uk2 )-a(T k ]^ k2 ) + (l 



aKa 2 kj8{kj - k$) (T h 



where, as previously, L, = \og(k 2 /kfy. 

Although this form of the equation does not look that much 
simpler at first sight than the original one, it has the advantage 
that it can be rewritten under the form of a Langevin equation 7 



The fluctuation term involves an awkward integration. By integration by 
part one can express it as a vertex applied to T but the integration over the 
internal variable z in |SJ is unknown. 
7 The complete hierarchy using |8} can also be rewritten as a Langevin equa- 
tion 1 1 2tl - Unfortunately, its expression is highly nontrivial. 
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FIG. 8: Speed of the wave for different values of the noise strength 
K(X 2 . The BK speed has been added for comparison. One sees that 
the speed decreases when K(X 2 increases. 
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FIG. 9: Numerical simulations of equation 1 141 . The dashed curves 
are different events corresponding to the same initial condition with 
different realisations of the noise term. Travelling waves are ob- 
served for each curve, together with dispersion. The black curve 
is the average amplitude. 



d Y T{L)=GL 



X(-d L )T(L) - T 2 (L) + Jia%T(L)l\(L,Y) 



(14) 

where T| is a Gaussian white noise satisfying the following 
commutation relations: 



{r\(L,Y)) 



0. 

4, 



(Ti(Li,yi)Ti(L2,y2)> = ^(L 1 -L 2 )8(Yi-Y 2 ). 



(15) 



The Langevin equation is a stochastic equation describing 
an event-by-event picture. Each realisation of the noise term 
corresponds to a particular evolution of the target and one can 
show that, once we average over those realisations using the 
correlations dl5> . the complete hierarchy for the evolution of 
(T k ) is recovered. It is interesting to notice that equation (1 1 4i 
is formally equivalent to the BK equation with an additional 
noise term. 

As for the case of the BK equation, one can restrict our- 
selves to the diffusive approximation. This leads to the 
stochastic FKPP (sFKPP) equation which amounts to add a 
noise term to the FKPP equation. It has found many appli- 
cations in various fields, e.g. in the description of reaction- 
diffusion systems. The noise term appears once we have to 
consider discreteness effects (e.g. for a reaction-diffusion sys- 
tem on a lattice, only an integer number of particles are ad- 
mitted per site). When the number of particles involved goes 
to infinity, the mean-field approximation is justified and the 
(non-stochastic) FKPP equation is valid. 

The basic effects of the additional noise term in the sFKPP 
equation are known, at least qualitatively, despite the fact that 
the fluctuations are only expected to have large consequences 
on the dilute tail of the wavefront, these modifications changes 
quite a lot the picture at saturation. In order to test the validity 
of those results in the case beyond the diffusive approxima- 
tion, i.e. with the full BFKL kernel, we have performed lll8ll 
numerical studies (see also lH^O of the QCD equation H4t . 
Moreover, many analytical results are only known in the limit 
where the noise strength Ka 2 is (irrealistically) small while 



we have concentrated our numerical work of physically ac- 
ceptable values. In the following paragraphs, we review the 
main effects of the fluctuation contribution and show that they 
are observed in the numerical analysis. 

First of all, for a single event, the evolved amplitude shows 
a travelling-wave pattern (up to small fluctuations in the far 
tail which, at least a this level, are irrelevant for the wave pat- 
tern at saturation). This means that each single realisation of 
the noise leads to geometric scaling, as it is the case for the 
BK equation. At this stage, the major difference comes from 
a decrease of the speed of the wave. Analytically, this can only 
be computed in the limit Ka 2 — > i.e. when the fluctuations 
are a small perturbation around the mean-field behaviour 8 . It 
has been found 120] that 



Vc 



aK 2 y c %"(y c ) 
21og 2 (a 2 K) 



Unfortunately, this expression gives only reliable results for 
extremely small values of Ka 2 — > such as 10~ 20 . This is 
clearly not sufficient for realistic situations. The numerical 
simulations we have performed clearly exhibits this decrease 
of the speed when the noise strength koc 2 increases (see figure 

The second noticeable effect of the noise term is to intro- 
duce dispersion between different events. Different events 
have the same shape but their position X(Y) (log(Q%(Y)) in 
physical units) fluctuates. At a give rapidity, one can compute 
the dispersion of those events. This diffusion process being 
very similar to a random walk, one expects that the dispersion 
of the events behaves like \/Y : 



(X 2 )y-(X)t 



D dis Y. 



' A recent analysis has also computed its strong-noise behaviour 12 m 
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FIG. 10: Dispersion (squared) of the position of the events as a func- 
tion of rapidity. As expected, it dispersion increases like yY but few 
dispersion is obtained in early stages of the evolution. 
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FIG. 11: Evolution of the averaged amplitude for different values 
of the noise strength, compared with the mean field result (dot- 
ted curve). The curves correspond, from left to right, to Y = 
0,10,20,30,40 and 50. 



The parametric dependence of the diffusion coefficient Ddiff 
has been computed numerically 1 20] and behaves 9 like 
|log- 3 (K0C?)|. 

Concerning the numerical simulations, the dispersion of the 
wavefront between different events and its increase with ra- 
pidity are manifest on figure [9] One can then extract the dis- 
persion as a function of rapidity. The resulting curve, shown 
on figure ^| for different values of the noise strength, calls 
for two remarks. First, the expected asymptotic behaviour 
(dispersion °c >/F) is observed and the diffusion coefficient 
increases with rapidity. However, for small values of the ra- 
pidity, we do not observe a significant dispersion. 

This dispersion of the events has an important physical 
consequences: although each single event displays geometric 
scaling, once we compute the average amplitude, dispersion 
induces geometric scaling violations 10 . This effect is best seen 



9 Recently, this coefficient has been computed analytically 22]. 
The asymptotic behaviour exhibits a new form of scaling, called diffusive 
scaling and its implications on DIS has been recently studied l23ll . 



on figure^2 wnere we nave compared the evolution with fluc- 
tuations to the BK results. While in the BK equation, a fixed 
travelling-wave pattern is formed, once fluctuations are take 
into account, we see a broadening of the average amplitude as 
rapidity increases. 

V. DISCUSSION AND PERSPECTIVES 

Let us now summarise the main points raised in these pro- 
ceedings. As this document is in itself some kind of a sum- 
mary, I shall only pick up the main points and refer to various 
references for detailed approaches. 

First, we have shown that it is possible to describe the evo- 
lution to high energy in pQCD by a hierarchy of equations 
(see (|8}). This gives the evolution of the (T) matrix (aver- 
aged over the target wave-function) and its higher-order corre- 
lations (r ). In this hierarchy, the evolution of (T k ) contains 
three types of contribution: 

1. the linear BFKL growth, proportional to d(T \. This 
is the usual high-energy contribution computed thirty 
years ago and corresponding to the exchange of k pQCD 
pomerons. It leads to a fast (exponential) increase of the 
scattering amplitude. 

2. the saturation corrections, proportional to a(r i+1 ). 
This negative term becomes important when the ampli- 
tude reaches unitarity and it allows for the constraint 
T(r,b) < 1 to be satisfied. 

3. the fluctuations term, proportional to ota^(r i_1 ). 
These fluctuations correspond to gluon-number fluctu- 
ations in the target. They play an important role when 
the amplitude is of order ocj or, equivalently, when the 
dipole density is of order one, i.e. in the dilute regime 
where one indeed expects fluctuation effects to appear. 

If one neglects the fluctuation contributions, we recover the 
Balitsky hierarchy (in its large-M- formulation) and, in the 
mean field approximation, the BK equation is obtained. This 
last equation, much simpler than the infinite hierarchy cap- 
tures many important points concerning the physics of satura- 
tion. We also have to point out that, as long as we do not take 
into account fluctuations, the evolution at all orders in 1 /N c 
exists under the form of the Balitsky/JIMWLK equation. We 
have only given a short introduction to the Colour Glass Con- 
densate (in which the JIMWLK equation is derived). Further 
information concerning the Balitsky (resp. CGC) approach 
can be found if reference ll24Tl (resp. ll25lo . 

Many additional things can be said concerning those dif- 
ferent contributions, e.g. concerning the equivalence between 
descriptions from the target and projectile point of view (see 
e.g. 12610 . or the attempts to describe the evolution beyond its 
large— A^ c - limit |27]. Again, the interested reader is forwarded 
to the corresponding references (and the references therein) 
for further information. 

In the second part of this overview, we have concentrated 
ourselves on the physical consequences arising from satura- 
tion and fluctuation effects. We have searched how those con- 
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tributions manifest themselves on the scattering amplitude and 
change the exponential behaviour obtained from the BFKL 
evolution. Again, the most important points can be sum- 
marised in two major steps with their respective physical con- 
sequences: 

1. Considering BFKL with saturation effects (mean-field 
picture), one can show that the BK equation lies in the 
same universality class as the FKPP equation. This im- 
plies that travelling waves are formed during the evolu- 
tion towards high energy. From the BFKL kernel one 
obtains the critical parameters which are the anomalous 
dimension observed in the large-^ 2 tail and the speed of 
the wave, equivalent to the exponent of the saturation 
scale. In physical terms, these travelling waves corre- 
sponds to geometric scaling. The experimental obser- 
vation at HERA can thus be seen as a consequence of 
saturation. It is extremely important to realise that ge- 
ometric scaling is a prediction from saturation physics 
which extends to large values of Q 2 , far beyond the sat- 
uration scale itself 11 . This geometric scaling has also 
been predicted at nonzero momentum transfer, a result 
which may be applied for example to p-meson produc- 
tion or DVCS. 

2. If one includes the effects of fluctuations, the evolu- 
tion becomes a Langevin equation 12 equivalent to the 
stochastic FKPP equation, including a noise term. For 
each realisation of this noise, we observe geometric 
scaling with a decrease of the speed w.r.t. the mean- 
field case. If we consider a bunch of events, we observe 
that, although they all have the same shape, their po- 
sition is diffused with diffusion going like y/Y. This 
dispersion of the saturation scale for a set of events im- 
plies that geometric scaling is violated for the averaged 
amplitude. Numerical studies indicate that geometric 
scaling is still valid at small rapidities. At higher ener- 
gies, a new (diffusive) scaling is predicted 12311 in which 
the amplitude scales like log(fc 2 /g 2 (T))/\/F. 

These two pictures, although they endow most of the phys- 

1 1 The geometric scaling window grows, in logarithmic units, like a/F beyond 
the saturation scale. 

12 If one consider the complete evolution equation with full phase-space de- 
pendence, one can also rewrite the hierarchy under the form of a Langevin 



ical effects, are far from being totally understood. Among the 
improvements, one can quote phenomenological applications 
like geometric ad diffusive scaling predictions for vector- 
meson productions, for diffraction 1 23] and for LHC physics. 
On more theoretical grounds most of the work needed to be 
done concerns the effects of fluctuations. Most of the results 
known so far are derived within the local approximation for 
the noise term and neglecting the impact parameter depen- 
dence. In addition, apart from numerical studies, analytical 
results can only be applied to irrealistically small values of 
CXj and a better understanding for more physical values is still 
lacking. Although the existing picture is expected to contain 
the qualitative effects, a more precise treatment is certainly an 
interesting problem to address. 

Within the framework presented here, the link between 
the QCD evolution equations and equations from statistical 
physics has led to a large number of interesting ideas and 
results. Even in statistical physics, many questions are yet 
opened, hence, extending the QCD picture including more so- 
phisticated treatments of the BFKL kernel and more precise 
forms of the noise term is certainly not a straightforward task. 
The questions of making clear predictions for phenomenol- 
ogy, of understanding the effects of fluctuations beyond the 
local-noise approximation and of including impact-parameter 
dependence are hence expected to give interesting work in the 
near future. Among particle physics, high-energy QCD is thus 
one of the most active fields. Because of those open questions 
and because of the requirement for an excellent knowledge 
of QCD at the LHC, we expect that high-energy QCD is still 
going to be active in the near future. 
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equation. It is also possible 28] to view the evolution s a reaction-diffusion 
process. 
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